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1.0  SUMMARY 


Aerospace  medicine  practitioners  track  the  epidemiology  of  conditions  that  limit  aircrew 
availability  and  work  toward  prevention  of  these  conditions.  These  prevention  efforts  should 
focus  on  those  conditions  that  are  the  primary  driver  of  aircrew  non-availability.  The  purpose  of 
this  study  was  to  reuse  available  datasets  to  conduct  an  analysis  of  potential  predictors  of  U.S. 

Air  Force  aircrew  non-availability  in  terms  of  being  in  “duties  not  to  include  flying”  (DNIF) 
status.  This  study  was  a  retrospective  cohort  analysis  of  U.S.  Air  Force  aircrew  on  active  duty 
during  the  period  from  2003-2012.  Predictor  variables  included  age,  Air  Force  Specialty  Code 
(AFSC),  clinic  location,  diagnosis,  gender,  and  pay  grade.  The  response  variable  was  DNIF 
duration.  Nonparametric  methods  were  used  for  the  exploratory  analysis  and  parametric  methods 
were  used  for  model  building  and  statistical  inference.  Significant  associations  were  observed 
between  age,  AFSC,  clinic,  and  primary  diagnosis  category  and  expected  DNIF  duration.  While 
controlling  for  specific  diagnoses,  increasing  age  was  positively  associated  with  expected  DNIF 
duration.  Six  AFSCs  were  associated  with  an  increased  expected  DNIF  duration;  however,  these 
AFSCs  were  not  significant  drivers  of  DNIF  duration  based  on  the  Pareto  analysis.  There  was 
observed  variability  in  expected  DNIF  duration  based  on  clinic,  with  six  clinics  identified  as 
significant  DNIF  drivers  after  controlling  for  other  demographic,  occupational,  and  health 
factors.  Forty  of  389  diagnostic  categories  were  included  and  significant  in  the  final  model. 

Based  on  Pareto  analysis,  25  of  these  primary  diagnosis  categories  were  significant  DNIF  drivers 
relative  to  the  other  diagnoses:  reproductive/pregnancy-related  conditions,  mental  health 
conditions,  fractures  and  degenerative  joint  conditions,  cardiopulmonary  conditions,  ocular 
conditions,  thyroid  disorders,  migraine  headaches,  enteritis  and  colitis,  hernias,  and  renal  calculi. 
Specific  demographic  (i.e.,  age),  occupational  (i.e.,  AFSC),  and  health  (i.e.,  clinic  location  and 
primary  diagnosis  category)  factors  were  identified  that  were  significantly  associated  with 
expected  DNIF  duration.  Subsequent  research  should  focus  on  the  application  of  primary, 
secondary,  and  tertiary  prevention  measures  to  ameliorate  the  potential  impact  of  these  DNIF 
drivers  where  possible. 

2.0  INTRODUCTION 

Practitioners  of  clinical  medicine  are  trained  to  prevent,  diagnose,  and  treat  conditions 
that  alter  a  patient’s  physiology  and  functional  state  in  a  normal  environment.  Practitioners  of 
aerospace  medicine  must  also  understand  the  interaction  of  a  patient’s  normal  or  abnormal 
physiology  and  functional  state  within  the  mission  environment  and  the  resulting  impact  on 
overall  flight  safety  and  performance.  Accordingly,  in  managing  acute  and  chronic  illnesses,  the 
aerospace  medicine  practitioner  has  the  additional  duty  of  rendering  an  aeromedical  disposition, 
that  is,  an  occupational  medicine  determination  whether  a  particular  aircrew  member  is  “fit  to 
fly.”  Prudent  aerospace  medicine  practitioners  also  track  the  epidemiology  of  conditions  that 
limit  aircrew  availability  and  work  toward  prevention  of  these  conditions  [1].  Given  ever-present 
resource  constraints,  not  the  least  of  which  is  aerospace  medicine  practitioner  time,  prevention 
efforts  should  focus  on  those  conditions  that  are  the  primary  driver  of  aircrew  non-availability. 
Unfortunately,  there  is  scant  published  literature  on  this  subject  to  inform  the  aerospace  medicine 
practitioner. 

The  purpose  of  this  study  was  to  reuse  available  datasets  to  conduct  an  exploratory 
analysis  of  potential  predictors  of  U.S.  Air  Force  (USAF)  aircrew  non-availability  in  terms  of 
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being  in  “duties  not  to  include  flying”  (DNIF)  status.  The  following  hypotheses  guided  this 
study: 

•  Hr.  Demographic  factors,  to  include  age  and  gender,  are  associated  with  duration  of 
DNIF  status. 

•  Hr.  Occupational  factors,  to  include  Air  Force  Specialty  Code  (AFSC),  service 
component,  and  pay  grade,  are  associated  with  duration  of  DNIF  status. 

•  Hr.  Health  factors,  in  terms  of  diagnoses  and  clinic,  are  associated  with  duration  of  DNIF 
status. 

3.0  METHODS 

3.1  Study  Design 

This  study  was  conducted  under  a  human-use  protocol  approved  by  the  71 1th  Human 
Performance  Wing  Institutional  Review  Board.  A  waiver  of  informed  consent  of  participants  was 
granted  due  to  the  impracticality  of  obtaining  written  consent  from  each  participant  in  the  study 
population.  This  study  was  a  retrospective  cohort  analysis  of  USAF  aircrew  on  active  duty 
during  the  period  from  2003-2012.  This  study  reused  a  dataset  created  for  a  study  analyzing  all 
outpatient  healthcare  encounters  occurring  in  any  of  the  USAF’s  Flight  and  Operational 
Medicine  Clinics  (FOMCs)  during  the  period  from  2003-2012  [2]  as  well  as  archival  data  on 
DNIF  events  extracted  from  the  Aeromedical  Services  Information  Management  System 
(ASIMS).  Inclusion  criteria  were  USAF  service  members  receiving  care  at  an  FOMC  with  at 
least  one  DNIF  episode.  Participants  were  excluded  if  they  had  missing  data  in  the  response 
variable,  DNIF,  or  in  the  personal  identifier. 

3.2  Data  and  Variables 

The  basic  unit  of  analysis  was  a  DNIF  episode.  The  duration  of  the  DNIF  episode  and  the 
associated  primary  diagnosis,  recorded  in  terms  of  International  Classification  of  Diseases,  Ninth 
Revision,  Clinical  Modification  (ICD-9-CM)  diagnosis  codes,  were  obtained  from  ASIMS. 
Participant  age  (continuous),  gender  (categorical  with  2  levels),  pay  grade  (categorical  with  16 
levels),  AFSC  (categorical  with  270  levels  [using  career  group,  career  field,  and  career  field 
subdivision  for  enlisted  personnel  and  career  group  and  functional  area  for  officers]),  service 
component  (categorical  with  3  levels),  and  FOMC  location  (categorical  with  77  levels)  for  each 
DNIF  episode  were  obtained  from  the  preexisting  study  database;  details  on  the  creation  of  this 
dataset  are  available  elsewhere  [2] . 

Diagnosis  codes  were  recoded  using  a  software  tool  developed  as  part  of  the  Healthcare 
Cost  and  Utilization  Project  (HCUP).  The  Clinical  Classification  Software  (CCS)  for  ICD-9-CM 
aids  analysts  to  collapse  diagnostic  data  from  over  14,000  diagnosis  codes  that  make  up  the 
ICD-9-CM  standardized  coding  system  into  clinically  meaningful  categories  [3].  The  367  tertiary 
level  classifications  were  used,  with  22  additional  levels  for  Department  of  Defense  (DoD) 
specific  categories,  such  as  “Medication  Education,”  “Armed  Forces  Health  Exam,”  and  “Travel 
Medication  Education,”  for  a  total  of  389  levels. 
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3.3  Statistical  Analysis 


The  original,  reused  dataset  [2]  comprised  90,331  distinct  participants.  A  total  of  7858 
participants  did  not  meet  the  study  inclusion  criteria  or  were  excluded  because  of  missing  data. 
The  final  study  population  comprised  389,976  DNIF  events  from  82,473  distinct  participants. 
The  study  dataset  was  randomly  partitioned  into  several  samples:  a  learning  sample  (235,919 
DNIF  events  from  50,000  distinct  participants)  for  exploratory  analysis  and  initial  variable 
selection,  a  training  sample  (70,150  DNIF  events  from  15,000  distinct  participants)  for  further 
variable  selection  using  a  marginal  longitudinal  model,  and  a  validation  sample  (71,938  DNIF 
events  from  15,000  distinct  participants)  for  statistical  inference  using  the  variables  selected  in 
the  first  two  steps.  A  fourth  remainder  sample  of  1 1,969  DNIF  events  from  2473  distinct 
participants  was  unused;  this  sample  was  held  back  in  case  further  data  exploration  was 
necessary.  Nonparametric  methods  were  used  for  the  exploratory  analysis  and  parametric 
methods  were  used  for  model  building  and  statistical  inference  given  the  greater  ease  of 
interpretation  of  the  latter  (e.g.,  standard  errors,  p-values,  etc.).  Separating  variable  selection  and 
model  building  ensured  that  the  reported  standard  errors  and  /;- values  were  valid. 

Tree-based  gradient  boosting  machine  (GBM)  [4]  modeling  was  used  for  exploratory 
analysis  on  the  learning  sample.  The  GBM  variable  importance  capability  was  used  to  select  the 
most  influential  predictors  to  include  in  the  parametric  analyses;  larger  variable  importance 
scores  suggested  greater  importance  in  terms  of  predicting  the  response.  Prior  to  analysis,  all 
high-level  categorical  variables  were  one-hot  encoded;  that  is,  a  separate  dummy  variable  was 
created  for  every  level  of  each  variable.  This  procedure  yielded  a  total  of  783  predictor  variables 
that  were  used  for  exploratory  analysis  of  the  GBM.  Variables  with  non-zero  importance  scores 
were  subsequently  included  in  the  parametric  analyses.  Since  the  study  objective  was  to  identify 
population-wide  predictors  of  DNIF  duration  rather  than  inference  on  individuals,  a  marginal 
model  rather  than  a  longitudinal  model  was  used  [5].  A  negative  binomial  model  with  a  log  link 
function  was  chosen  because  the  response  variable  was  a  count  variable  with  dissimilar  mean 
and  variance  (thus  making  a  Poisson  model  a  suboptimal  choice).  Predictor  variables  included 
age,  gender,  pay  grade,  AFSC,  service  component,  FOMC  location,  and  diagnosis.  Participant 
was  a  random  repeated  measure  in  the  model,  and  a  compound  symmetry  (exchangeable) 
covariance  structure  was  assumed. 

R  version  3.3.2  [6]  was  used  for  data  preparation  and  calculation  of  summary  statistics. 
The  R  gbm  package,  version  2.1.1,  was  used  to  accomplish  the  GBM  modeling.  SAS  version  9.4 
(SAS  Institute,  Cary,  NC)  was  used  to  create  the  sample  datasets  (PROC  SURVEYSELECT). 
SAS  (PROC  GENMOD)  was  used  to  fit  the  marginal  longitudinal  model  on  the  training  sample 
and  estimate  the  model  of  the  validation  sample.  Statistical  significance  was  defined  as 
p  =  0.0001. 
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4.0  RESULTS 


4.1  Descriptive  Statistics 

Table  1  provides  descriptive  statistics  for  the  measured  variables  for  the  final  study 
population.  Clinic  location  is  not  displayed,  as  only  clinic  pseudoidentifiers  were  provided  to 
preserve  data  de-identification.  With  the  exception  of  primary  diagnosis  category,  summary 
statistics  were  computed  on  the  basis  of  the  population  of  unique  participants  using  a  randomly 
selected  DNIF  event  during  the  first  year  of  observation  to  establish  a  measurement  for  each 
variable.  In  contrast,  summary  statistics  for  primary  diagnosis  category  were  computed  based  on 
the  population  of  DNIF  events. 

4.2  Hypothesis  Testing 

In  the  final  validation  sample,  the  predicted  results  of  the  marginal  longitudinal  model  for 
DNIF  duration  had  a  correlation  (r)  of  0.45  with  the  actual  number  of  DNIF  days.  Out  of  the  783 
predictor  variables  used  in  the  GBM  model  fitted  on  the  learning  data,  339  variables  had  a  non¬ 
zero  relative  influence  and  were  included  in  the  parametric  analyses.  Of  these  predictor 
variables,  84  variables  exclusive  of  the  intercept  had  statistically  significant  associations  at 
p  <  0.0001  with  DNIF  duration  when  the  initial  negative  binomial  model  was  fitted  using  the 
training  data.  Fifty-two  variables,  not  including  the  intercept,  had  statistically  significant 
associations  at  alpha  p  <  0.0001  with  DNIF  duration  when  the  final  negative  binomial  model  was 
fitted  on  the  validation  dataset  (Table  2). 

Based  on  the  model  results,  we  partially  accept  hypothesis  1  that  demographic  factors  are 
associated  with  the  duration  of  DNIF  status.  There  was  a  significant  association  with  age  and 
duration  of  DNIF  status,  while  gender  was  not  a  predictor  variable  selected  for  inclusion  in  the 
model.  We  also  partially  accept  hypothesis  2  that  occupational  factors  are  associated  with  the 
duration  of  DNIF  status.  Of  the  occupational  factors  considered,  only  AFSC  was  selected  for 
inclusion  in  the  model,  and  then  only  6  of  the  potential  270  levels  of  this  variable  were  included 
and  significant  in  the  final  model.  We  accept  hypothesis  3  that  health  factors  are  associated  with 
the  duration  of  DNIF  status.  Forty  of  389  diagnostic  categories  were  included  and  significant  in 
the  final  model,  while  7  out  of  77  potential  clinic  locations  were  included.  Figure  1  provides  a 
Pareto  display  of  the  primary  diagnosis  categories  that  were  significantly  associated  with 
expected  days  DNIF.  Six  clinics  and  25  diagnosis  categories  were  the  primary  drivers  of  DNIF 
duration  based  on  observed  effect  size. 
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Table  1.  Descriptive  Statistics  for  the  Study  Population 


Variable 

Descriptive  Statistic 

N 

82,473 

DNIF  events: 

N 

389,976 

Duration,  days,  median  (IQR) 

7.00(13.58) 

Age,  yr,  median  (IQR) 

27  (12) 

Gender  (ref  male),  no.  (%) 

72,  834  (88.31) 

Service  component,  no.  (%): 

Active  duty 

75,379  (91.40) 

Reserve 

3909  (4.74) 

National  Guard 

3185  (3.86) 

AFSC,  no.  (%): 

Officer: 

92TX  Pilot  trainee 

7333  (20.59) 

11FX  Fighter  pilot 

3351  (9.41) 

13SX  Space  and  missile 

2824  (7.93) 

1 1MX  Mobility  pilot 

2697  (7.57) 

1 1  AX  Airlift  pilot 

2068  (5.82) 

11KX  Trainer  pilot 

1747  (4.91) 

13BX  Air  battle  manager 

1570  (4.41) 

11TX  Tanker  pilot 

947  (2.67) 

11RX  Reconnaissance/surveillance/electronic  warfare  pilot 

784  (2.20) 

12RX  Reconnaissance/surveillance/electronic  warfare  combat  systems  officer 

783  (2.20) 

12BX  Bomber  combat  systems  operator 

713(2.00) 

1 1BX  Bomber  pilot 

629(1.77) 

11SX  Special  operations  pilot 

563  (1.58) 

12FX  Fighter  combat  systems  officer 

561  (1.57) 

62EX  Developmental  engineer 

557(1.56) 

46FX  Flight  nurse 

486(1.36) 

1 1FIX  Flelicopter  pilot 

481  (1.35) 

92SX  Student  officer  authorization 

476(1.34) 

Other3 

7040  (19.77) 

Enlisted: 

1C1XX  Air  traffic  control 

4773  (13.06) 

1A2XX  Aircraft  loadmaster 

4233  (11.58) 

9T0XX  Basic  enlisted  airman 

2764  (7.56) 

1A1XX  Flight  engineer 

2344  (6.41) 

1A8XX  Airborne  cryptologic  linguist 

2303  (6.30) 

1A3XX  Airborne  mission  system 

1718  (4.70) 

1C6XX  Space  systems  operations 

1469  (4.02) 

1A0XX  In-flight  refueling 

1447  (3.96) 

4N0XX  Aerospace  medical  service 

1313  (3.59) 

1A4XX  Airborne  operations 

1259  (3.44) 

1C5XX  Command  and  control  battle  management  operations 

1118(3.06) 

2A5XX  Aerospace  maintenance 

818  (2.24) 

3P0XX  Security  forces 

802  (2.19) 

1C4XX  Tactical  air  control  party 

749  (2.05) 

1C2XX  Combat  control 

730  (2.00) 

1T2XX  Pararescue 

729(1.99) 

1N1XX  Geospatial  intelligence 

569(1.56) 
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Table  1.  Descriptive  Statistics  for  the  Study  Population  (concluded) 


Variable 

Descriptive  Statistic 

Enlisted  AFSC  (continued): 

1A7XX  Aerial  gunner 

513  (1.4) 

Other3 

6906(18.91) 

Missing 

10,306(12.50) 

Primary  diagnosis  category, b  no.  (%): 

Diseases  of  the  respiratory  system 

104,637  (26.83) 

DoD  specific:  education  or  counseling 

48,117  (12.34) 

Diseases  of  the  digestive  system 

31,177  (7.99) 

Diseases  of  the  nervous  system  and  sense  organs 

30,625  (7.85) 

Symptoms;  signs,  ill-defined  conditions  and  factors  influencing  health  status 

26,360  (6.76) 

Diseases  of  the  musculoskeletal  system  and  connective  tissue 

24,521  (6.29) 

Injury  and  poisoning 

22,404  (5.74) 

Diseases  of  the  skin  and  subcutaneous  tissue 

5529(1.42) 

Infectious  and  parasitic  diseases 

5425  (1.39) 

Residual  codes,  unclassified,  all  E  codes 

5228(1.34) 

Complications  of  pregnancy;  childbirth;  and  the  puerperium 

4917(1.26) 

Diseases  of  the  genitourinary  system 

4899(1.26) 

Diseases  of  the  circulatory  system 

4277  (1.10) 

Mental  illness 

3494  (0.90) 

Endocrine;  nutritional;  and  metabolic  diseases  immunity  disorders 

3129  (0.80) 

DoD  specific  exams 

1228  (0.31) 

Neoplasms 

768  (0.20) 

Congenital  anomalies 

462  (0.12) 

Diseases  of  the  blood  and  blood-forming  organs 

174  (0.04) 

Certain  conditions  originating  in  the  perinatal  period 

83  (0.02) 

Other  DoD  specific  diagnoses 

16  (<0.01 ) 

DoD  specific:  traumatic  brain  injury 

5  (<0.01) 

Missing 

62,501  (16.03) 

Pay  grade,  no.  (%): 

Officer 

01 

16,433  (36.29) 

02 

5180(11.44) 

03 

11,057  (24.42) 

04 

6045  (13.35) 

05 

5178(11.44) 

06+ 

1388  (3.07) 

Enlisted 

El 

1688  (4.54) 

E2 

2366  (6.36) 

E3 

10,809  (29.06) 

E4 

5152  (13.85) 

E5 

8319  (22.37) 

E6 

4935  (13.27) 

E7 

2979  (8.01) 

E8 

753  (2.02) 

E9 

191  (0.51) 

Note:  IQR  =  interquartile  range. 

“Only  AFSCs  comprising  80%  of  participants  shown  for  brevity. 
bHCUP-CCS  secondary  level  diagnosis  categories  shown  for  brevity. 
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Table  2.  Negative  Binomial  Regression  Model  Results  for  DNIF  Duration 


Variable 

B 

SE  (B) 

Expected 
Days  DNIF 

p -value 

Intercept 

2.7409 

0.0459 

15.50 

<0.0001 

Age 

0.0219 

0.0013 

15.84 

<0.0001 

Gender  (ref  =  male) 

0.1019 

0.0280 

17.16 

0.0003 

Clinic  location: 

CLID025606592 

0.0950 

0.0564 

17.05 

0.0925 

CLID068902320 

0.0866 

0.0648 

16.90 

0.1810 

CLID093047257 

0.2947 

0.0742 

20.81 

<0.0001 

CLID 106868943 

0.2582 

0.0673 

20.07 

0.0001 

CLID 1096 16999 

0.2022 

0.0616 

18.97 

0.0010 

CLID  110435376 

0.6035 

0.1325 

28.34 

<0.0001 

CLID1 11016652 

-0.5729 

0.0517 

8.74 

<0.0001 

CLID  125682959 

0.4580 

0.0459 

24.51 

<0.0001 

CLID  14785 1280 

0.4382 

0.1614 

24.03 

0.0066 

CLID254322565 

0.3939 

0.0644 

22.98 

<0.0001 

CLID322301264 

0.1306 

0.0682 

17.66 

0.0555 

CLID381735261 

0.8947 

0.1019 

37.92 

<0.0001 

Primary  diagnosis  category: 

Acute  bronchitis 

-0.7744 

0.0811 

7.15 

<0.0001 

Administrative/social  admission 

-0.3520 

0.0919 

10.90 

0.0001 

Allergic  reactions 

-0.3797 

0.1013 

10.60 

0.0002 

Bipolar  disorders 

1.3169 

0.1711 

57.85 

<0.0001 

Blindness  and  vision  defects 

0.3070 

0.0505 

21.07 

<0.0001 

Calculus  of  urinary  tract 

0.5757 

0.1134 

27.57 

<0.0001 

Cardiac  dysrhythmias 

1.0505 

0.1189 

44.32 

<0.0001 

Cataract 

1.0928 

0.2491 

46.23 

<0.0001 

Cellulitis  and  abscess 

-0.8052 

0.0789 

6.93 

<0.0001 

Codes  related  to  mental  health  disorders 

1.0261 

0.1731 

43.25 

<0.0001 

Coronary  atherosclerosis  and  other  heart  disease 

0.9879 

0.2387 

41.63 

<0.0001 

Depressive  disorders 

1.3929 

0.1622 

62.41 

<0.0001 

DoD  specific:  medication  education 

-0.8379 

0.0274 

6.71 

<0.0001 

Ectopic  pregnancy 

1.6414 

0.0956 

80.02 

<0.0001 

Encephalitis,  except  that  caused  by  TB  or  STD 

-0.8920 

0.3284 

6.35 

0.0066 

Endometriosis 

0.5373 

0.3597 

26.53 

0.1352 

Essential  hypertension 

0.5179 

0.0769 

26.02 

<0.0001 

Fracture  of  lower  limb 

0.9274 

0.0780 

39.19 

<0.0001 

Fracture  of  upper  limb 

0.7325 

0.0566 

32.25 

<0.0001 

Gastritis  and  duodenitis 

-0.8559 

0.1386 

6.59 

<0.0001 

Glaucoma 

0.5312 

0.1635 

26.37 

0.0012 

Heart  valve  disorders 

0.7483 

0.3682 

32.76 

0.0421 

Influenza 

-0.9445 

0.1113 

6.03 

<0.0001 

Inguinal  hernia 

0.4951 

0.0956 

25.43 

<0.0001 

Intervertebral  disc  disorders 

1.2521 

0.0769 

54.22 

<0.0001 

Migraine 

1.1869 

0.1766 

50.80 

<0.0001 

Nausea  and  vomiting 

-1.1883 

0.0951 

4.72 

<0.0001 

Non-Hodgkins  lymphoma 

0.2085 

0.2727 

19.09 

0.4444 

Other  abdominal  hernia 

0.4652 

0.1009 

24.68 

<0.0001 

Other  aftercare 

0.1566 

0.0426 

18.13 

0.0002 

Other  and  ill-defined  heart  disease 

1.5975 

0.5827 

76.58 

0.0061 
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Table  2.  Negative  Binomial  Regression  Model  Results  for  DNIF  Duration  (concluded) 


Variable 

B 

SE(B) 

Expected 
Days  DNIF 

p -value 

Primary  diagnosis  category  (continued): 

Other  and  unspecified  gastrointestinal  disorders 

-0.7095 

0.1148 

7.62 

<0.0001 

Other  and  unspecified  asthma 

1.3220 

0.4394 

58.14 

0.0026 

Other  chronic  pulmonary  disease 

-0.9553 

0.0708 

5.96 

<0.0001 

Other  complications  of  pregnancy 

1.3600 

0.1821 

60.39 

<0.0001 

Other  fractures 

0.7649 

0.0961 

33.31 

<0.0001 

Other  mycoses 

-0.7298 

0.1116 

7.47 

<0.0001 

Other  non-traumatic  joint  disorders 

0.3820 

0.0557 

22.71 

<0.0001 

Other  thyroid  disorders 

1.3024 

0.1407 

57.01 

<0.0001 

Other  upper  respiratory  infections 

-1.0044 

0.0239 

5.68 

<0.0001 

Other  viral  infections 

-0.9378 

0.0642 

6.07 

<0.0001 

Otitis  media  and  related  conditions 

-0.6187 

0.0689 

8.35 

<0.0001 

Outcome  of  delivery  (V  codes) 

1.7020 

0.0575 

85.02 

<0.0001 

Peri-;  endo-;  &  myocarditis;  cardiomyopathy  (except  that 

0.4334 

0.3284 

23.91 

0.1869 

caused  by  TB  or  STD) 

Phlebitis;  thrombophlebitis  and  thromboembolism 

1.5985 

0.3762 

76.66 

<0.0001 

Pneumonia  (except  that  caused  by  TB  or  STD) 

-0.4967 

0.1524 

9.43 

0.0011 

Pulmonary  heart  disease 

1.2864 

0.3844 

56.11 

0.0008 

Regional  enteritis  &  ulcerative  colitis 

1.3004 

0.3220 

56.90 

<0.0001 

Residual  codes;  unclassified;  all  E  codes 

0.2626 

0.0660 

20.16 

<0.0001 

Retinal  detachments;  defects;  vascular  occlusion;  and 

0.8921 

0.1102 

37.83 

<0.0001 

retinopathy 

Spondylosis  and  allied  disorders 

1.0224 

0.3723 

43.09 

0.0060 

Sterilization 

-0.7122 

0.0604 

7.60 

<0.0001 

Urinary  tract  infections 

-0.9415 

0.1173 

6.05 

<0.0001 

AFSC: 

1 1EX  Experimental  test  pilot 

-0.2035 

0.2479 

12.65 

0.4116 

1 1FX  Fighter  pilot 

-0.3321 

0.0444 

11.12 

<0.0001 

1 1KX  Trainer  pilot 

-0.2595 

0.0498 

11.96 

<0.0001 

1 1MX  Mobility  pilot 

-0.2355 

0.0302 

12.25 

<0.0001 

1 1RX  Reconnaissance/surveillance/electronic  warfare  pilot 

-0.2932 

0.0611 

11.56 

<0.0001 

1C  1X1  Air  traffic  control 

-0.3921 

0.0379 

10.47 

<0.0001 

1C2X0  Combat  control 

0.4644 

0.3949 

24.66 

0.2396 

1C3X1  Command  post 

-0.0124 

0.4911 

15.31 

0.9798 

1N3X4  Cryptologic  language  analyst 

0.7745 

0.3597 

33.63 

0.0313 

21RX  Logistics  readiness 

0.2455 

0.3964 

19.81 

0.5358 

3C0X1  Communication-computer  systems 

-0.2341 

0.3440 

12.27 

0.4961 

3E5X1  Engineering 

-0.5207 

0.2859 

9.21 

0.0686 

44AX  Chief,  hospital/clinic  services 

0.6519 

0.2641 

29.75 

0.0136 

48AX  Aerospace  medicine  specialist 

-0.1738 

0.1555 

13.03 

0.2637 

4E0X1  Public  health 

-1.1206 

0.1636 

5.05 

<0.0001 

83RX  Recruiting  service 

0.0118 

0.2995 

15.68 

0.9685 

91WX  Wing  commander 

-0.3013 

0.1595 

11.47 

0.0590 

Note:  CLID  =  clinic  identifier;  SE  =  standard  error;  STD  =  sexually  transmitted  disease;  TB  =  tuberculosis. 
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Expected  Days  DNIF 

—  K)U»4^c/iOv-^JOO^O 
oooooooooo 


Outcome  of  delivery  (V  codes) 
Ectopic  pregnancy 
Phlebitis/thrombophlebitis/thromboembolism 
Depressive  disorders 
Other  complications  of  pregnancy 
Bipolar  disorders 
Other  thyroid  disorders 
Regional  enteritis  &  ulcerative  colitis 
Intervertebral  disc  disorders 
Migraine 
Cataract 
Cardiac  dysrhythmias 
Codes  related  to  mental  health  disorders 
Coronary  atherosclerosis  and  other  heart  disease 
Fracture  of  lower  limb 
CLID381735261 
Retinal  detachments;  defects;  vascular  occlusion; 

Other  fractures 
Fracture  of  upper  limb 
CLID1 10435376 
Calculus  of  urinary  tract 
Essential  hypertension 
Inguinal  hernia 
Other  abdominal  hernia 
CUD  125682959 
CLID254322565 
Other  non-traumatic  joint  disorders 
Blindness  &  vision  defects 
CLID093047257 
Residual  codes;  unclassified;  all  E  codes 
CLID 106868943 
Age 

1 1MX  Mobility  pilot 
11KX  Trainer  pilot 
1 1RX  Reconnaissance/surveillance/electronic 
11FX  Fighter  pilot 
Administrative/social  admission 
1C  1 X 1  Air  traffic  control 
CLID11 1016652 
Otitis  media  &  related  conditions 
Other  and  unspecified  gastrointestinal  disorders 

Sterilization 
Other  mycoses 
Acute  bronchitis 
Cellulitis  &  abscess 
DOD  specific:  medication  education 
Gastritis  and  duodenitis 
Other  viral  infections 
Urinary  tract  infections 
Influenza 

Other  chronic  pulmonary  disease 
Other  upper  respiratory  infections 
4E0X1  Public  health 
Nausea  &  vomiting 


Figure  1.  Pareto  display  of  the  significant  predictors  of  expected  DNIF  duration. 
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5.0  DISCUSSION 


To  the  best  of  the  authors’  knowledge,  this  study  is  the  first  attempt  to  systematically 
explore  potential  predictors  of  US  AF  aircrew  non-availability  in  terms  of  being  in  DNIF  status 
over  a  10-year  period.  Significant  associations  were  observed  between  age,  AFSC,  clinic,  and 
primary  diagnosis  category  and  expected  DNIF  duration.  While  controlling  for  specific 
diagnoses,  increasing  age  was  positively  associated  with  expected  DNIF  duration.  Six  AFSCs 
were  associated  with  an  increased  expected  DNIF  duration;  however,  these  AFSCs  were  not 
significant  drivers  of  DNIF  duration  based  on  the  Pareto  analysis.  There  was  observed  variability 
in  expected  DNIF  duration  based  on  clinic,  with  six  clinics  identified  as  significant  DNIF  drivers 
after  controlling  for  other  demographic,  occupational,  and  health  factors.  As  anticipated,  multiple 
primary  diagnosis  categories  were  associated  with  increased  expected  DNIF  duration.  Based  on 
Pareto  analysis,  25  of  these  diagnosis  categories  appeared  to  be  significant  DNIF  drivers  relative 
to  the  other  diagnoses:  reproductive/pregnancy-related  conditions,  mental  health  conditions, 
fractures  and  degenerative  joint  conditions,  cardiopulmonary  conditions,  ocular  conditions, 
thyroid  disorders,  migraine  headaches,  enteritis  and  colitis,  hernias,  and  renal  calculi.  Of  note, 
gender  was  not  associated  with  expected  DNIF  duration  after  controlling  for  diagnoses. 

Given  this  analysis,  the  next  step  is  to  evaluate  those  conditions  found  to  be  significant 
DNIF  drivers  and  identify  opportunities  for  primary,  secondary,  and  tertiary  prevention  [7]. 

Since  infectious  diseases  were  not  among  the  DNIF  drivers,  traditional  primary  prevention 
measures  focusing  on  vaccination  are  of  limited  utility.  Instead,  primary  prevention  should  focus 
on  those  conditions  caused  by  injuries  and/or  toxic  exposures  resulting  from  modifiable 
environmental  exposures.  Secondary  prevention  should  focus  on  screening,  either  for  specific 
conditions  or  antecedent,  modifiable  risk  factors  for  those  conditions  (e.g.,  hypertension  and 
coronary  atherosclerotic  disease).  Routine  screening  is  already  accomplished  as  part  of  the 
mandated,  annual  periodic  health  assessment.  Subsequent  research,  however,  is  needed  to 
correlate  current  screening  tools  with  observed  DNIF  drivers.  Finally,  tertiary  prevention 
activities  should  seek  to  minimize  expected  DNIF  duration  after  a  condition  occurs  by 
optimizing  treatment  selection  and  delivery  throughout  the  care  cycle  for  the  condition. 

In  conclusion,  specific  demographic  (i.e.,  age),  occupational  (i.e.,  AFSC),  and  health  (i.e., 
clinic  location  and  primary  diagnosis  category)  factors  were  identified  that  were  significantly 
associated  with  expected  DNIF  duration.  Subsequent  research  should  focus  on  the  application  of 
primary,  secondary,  and  tertiary  prevention  measures  to  ameliorate  the  potential  impact  of  these 
DNIF  drivers  where  possible. 
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LIST  OF  ABBREVIATIONS  AND  ACRONYMS 


AFSC  Air  Force  Specialty  Code 

ASIMS  Aeromedical  Services  Information  Management  System 

CCS  Clinical  Classification  Software 

CLID  clinic  identifier 

DNIF  duties  not  to  include  flying 

DoD  Department  of  Defense 

FOMC  Flight  and  Operational  Medicine  Clinic 

GBM  gradient  boosting  machine 

HCUP  Healthcare  Cost  and  Utilization  Project 

ICD-9-CM  International  Classification  of  Diseases,  Ninth  Revision,  Clinical  Modification 

IQR  interquartile  range 

SE  standard  error 

STD  sexually  transmitted  disease 

TB  tuberculosis 


USAF  U.S.  Air  Force 
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